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ABSTRACT 

Organic material found in meteorites and interplanetary dust particles is enriched 
in D and 15 N. This is consistent with the idea that the functional groups carrying these 
isotopic anomalies, nitriles and amines, were formed by ion-molecule chemistry in the 
protosolar nebula. Theoretical models of interstellar fractionation at low temperatures 
predict large enrichments in both D and 15 N and can account for the largest isotopic 
enrichments measured in carbonaceous meteorites. However, more recent measure- 
ments have shown that, in some primitive samples, a large 15 N enrichment does not 
correlate with one in D, and that some D-enriched primitive material displays little, 
if any, 15 N enrichment. By considering the spin-state dependence in ion-molecule re- 
actions involving the ortho and para forms of H 2 , we show that ammonia and related 
molecules can exhibit such a wide range of fractionation for both 15 N and D in dense 
cloud cores. We also show that while the nitriles, HCN and HNC, contain the greatest 
15 N enrichment, this is not expected to correlate with extreme D enrichment. These 
calculations therefore support the view that Solar System 15 N and D isotopic anoma- 
lies have an interstellar heritage. We also compare our results to existing astronomical 
observations and briefly discuss future tests of this model. 

Subject headings: Astrochemistry — ISM: molecules — ISM: clouds — Meteorites, 
meteors, meteoroids — Molecular processes — Planets and satellites: formation 



Introduction 



The isotopic enrichments measured in t he primitive organic matter found in comets, interplan - 
etary dust particles (IDPs) and meteorites (|Remusat et al.l 12006c Floss et al.l 120061 : lAleonl l2010h . 
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probably had their chemical origin in a low-temperature environment that faciliated ion-molecule 
isotopi c exchange reactions, such as the dense interstellar medium (ISM) or the outer protosolar 
nebula (|Mumma & CharnlewOl ll : [Messenger et al.ll2003l) . 



The enhanced D/H ratios seen in some primitive mater ial (|Messengerl l2000l : iRobertl 12003 



Messenger et al.ll2003l : Keller et al.ll2004l ; lHartogh et al.ll201 II) almost certainly originated in low- 
temperature enviro nments. Deuteriu m enhancement is initiated by the exothermic ion-molecule 
exchange reaction (|Millar et al.ll 19891) going to the right 



HD ^ H 2 D~ 



H, 



170 K 



(1) 



In dense cores where CO and other heavy molecules have condensed as ices on to dust grains 
(i.e. become depleted) , H 2 D + molecules can continue to react with HD, producing HD^ and 



D3 (|Roberts et al.l 120031) . This produces very high levels of deuteration in the remaining gaseous 
species as well as, th rough the resulting high atomic D/H ratios, in molecular ice mantles (e.g. 
Ceccarelli et al.1 12007). A similar D fractionation chemistry can also occur in the cold regions of 
protoplanetary disks (IWillac yll2Q07h . 



In both Jupiter Family and Port Cloud comets, N/ N ratios are measu red (~ 130 — 170) for 



both CN and HCN (IHutsemekers et al.ll2008l : iBockelee-Morvan et all 120081) that are signi ficantly 



enhan ced relative to both the terrestrial and protosolar values (272 and 440 respectively, iMarty 
20121) . Primitive materials can exhibit large bulk 15 N enhancements ( 14 N/ 15 N ~ 180 - 300), 



which exhibit the greatest D and 15 N enhancements known ( 14 N/ 15 N < 100) (IMessenger 


2000; 


Aleon et al. 


2003; 


Floss et al. 


2006; 


Busemann et al. 


2006; 


Briani et al. 


2009; 


Aleon 


2010). 



Theoretical models of nitrogen isotopic chemistry in molecular cloud cores, where gaseous 
molecules are being depleted by accretion on to d ust grains, ca n reproduce the range of 14 N/ 15 N 



Aleonl |2010|) throu gh exchange reactions be 



ratios found in primitive matter (~ 50 — 280, 

tween 15 N atoms and molec ular ions such as 14 N 2 H + and HC 14 NH + (|Charnley & Rodgersl 12002 ; 
Rodgers & Charnlevll2008allbl) . These models predict that N 2 , NH 3 , HCN, CN and HN C molecules 
should possess the greatest 15 N enrichments. There are two 15 N-fractionating pathways (|Rodgers & Charnley 
2008bl) : a slow one to N 2 and ammonia (~ 10 6 years) and more rapid one to HCN and other nitriles 
(~ 10 5 years). 

Recent, albeit sparse, astronomical observations of starless cores indicate that interstellar 
NH 3 and N 2 are not as enriched as p redicted for the ISM or as mea sured in the Solar Sy stem: 
14 NH 3 / 15 NH 3 « 334 Lis et alJboioR and 14 N 2 H+/ 14 N 15 NH+ « 446 iBizzocchi et alJboioh . 



However, the most pressing problem for a direct ISM-Solar System isotopic connection con- 
cerns the fact that, while 15 N and D meteoritic hotspots do seem to correlate in some samples (e.g. 
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Aleon et al.ll2003|) . they clearly do not in o thers (e.g. LBusemann et al J 120061: iGourier et al.ll20 08) 



Yet, as has been n oted by several authors ([Alexander et al J 12008c iBriani et al .1120091 : iBonal et al 



2010l : lAleonl |2010() . the low-temperature interstellar environments most conducive to producing 
large 15 N enrichments, should also produce concomitantly enormous molecular D/H ratios (e.g. 
in ammonia or HCN). Thus, one would expect 15 N and D hotspots to always be spatially corre- 
lated, contrary to what is seen. This presents a serious challenge for ion-molecule fractionation 
chemistry. 

In cold interstellar envi ronments, nuclear s pin states can be an important factor in determining 
chemical reaction rates (e.g. Flower et al. The major interstellar collision partner, molecular 

hydrogen, has a 170.5 K difference in zero point energy between the para (antiparallel spins) and 
the ortho (aligned nuclear spins) forms; this allows some endoergic reactions to pro ceed, even at 



low te mperatures, if a substantial fraction of the H 2 is in the higher energy form, o-H 2 . Pagani et al 



(|201ll) demonstrated that a high H 2 ortho-to-para ratio (OPR) therefore inhibits ('poisons') the 
production of H 2 D + in reaction (1), suppressing the overall deuteration. 

The fractionation of 15 N may also be influenced by the OPR of H 2 . Ammonia formation is 
initiated by the production of N + from N 2 by He + , which then reacts in 

N + + H 2 — ► NH+ + H, (2) 

followed by a sequence of ion-molecule reactions with H 2 through to NH t and a final electron 
dissociative recombination. Reaction |2] has an activation energy of <200 K (|Gerrichll 19931) which, 
at low temperatures, can be overcome by the internal energy of o-H 2 . Thus, a mmonia formatio n 
and the related 15 N fractionation in dark clo uds are dependent on the H 2 OPR (|Le Bourlotlll99ll) . 
a distinction not made in earlier models (e.g. iRodgers & Charnlevll2008bl). A recent re-assessment 
of the OPR dependence in the original experimental data by iDislaire et al.l (|2012l) indicates that 
previous work has overestimated the low-temperature rate coefficient for reaction [2l involving o- 
H 2 , by almost three orders of magnitude (cf. Le Bourlot 1991). Thus, the precise abundance of 
o-H 2 could play a pivotal role in producing a diverse range of D- 15 N fractionation in interstellar 
precursor molecules. 

In this Letter we quantify the effect the H 2 OPR has on interstellar 15 N fractionation, and 
demonstrate that this can account for the isotopic anomalies observed in primitive Solar System 
materials. 



2. Model 



W e consider the chemical evolution of the central regions of a static prestellar core (|Charnley & Rodgers 
2002|) with a gas density n(H 2 )=10 6 cm -3 , a temperature of 10 K and a visual extinction of Ay>10 
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magnitudes. Cos mic ray ionization oc c urs, at a rate of (=3 x 10 17 s 1 . The chemical model used 



here is based on Rodgers & Charnley ( 2008a b_U and Fig. Q] shows the most important reactions 
for 15 N f ractionation. Helium , carbo n, oxygen and nitrogen are included at elemental abundances 
given by ISavage & Sembachj (I1996I) . We are interested in cores that are just about to form pro- 
tostars and so we assume that all carbon is initially bound up in CO, with the remaining oxygen 
in atomic form, and the nitrogen is partly atomic with 50% in molecular and atomic form respec- 
tively. The elemental 14 N/ 15 N ratio is assumed to be 440, with the same fraction of 15 N in atomic 
form initially, 15 N/ 15 N 14 N=1. 



The gas-grain interaction is modelled according to ICharnleyl (|1997|) : all neutral gas-phase 
species, except H 2 , He, N, and N 2 , stick and freeze out onto grains upon collision. No grain 
surface chemistry or desorption is considered, except for sticking of H atoms (with an efficiency of 
0.6) followed by reaction to produce H 2 molecules which are ejected upon formation (see Sect |2.1|) . 



2.1. Inclusion of spin-state reactions 



The reaction network has been expanded to contain ortho and para forms of H 2 , H^ , and 
H^", and couples 279 gas-phase species by 4420 reactions. Both spin types have been consid- 
ered eq ually likely to react with a given species, with the exception of reaction ©, where new 
rates of bislaire etaH d2012T> for ortho and para H 2 , k^ o =4.20x 10~ 10 (T/300)~ ai7 exp(-44.5/T) 
cm 3 s _1 and h^ p =S.35x 10~ 10 exp(— 168. 5/T) cm 3 s" 1 , respectively, have been adopted. Other 
temperature-dependent H 2 reactions have energy barriers that are too high for the ground- state 
energy difference between the spin types to have a significant effect at 10 K. 

Initially the H 2 OPR is s et to the high tempera ture statistical equilibrium value of 3. Measure- 
ments in the diffuse ISM by ICrabtree et al.1 (|201ll) are consistent with an H 2 OPR around 1 ; we 
find adopting this lower value has little effect on our results. The reaction 



H + + o-H 2 



H" 



v- 



-H, 



(3) 



is exothermic by 170.5 K and so o-H 2 is efficiently converted into p-H 2 at low temperatures. We 
assume that the H 2 formed and ejected from grains has an OPR of 3; this ensures a steady-state 
OPR is reached. Reaction rates for conversion between ortho and para sp ecies, as well as for thei r 
interaction with electrons and cosmic rays, have been implemented from IWalmsley et al.l (|2004l) . 
Branching ratios fo r o/p-Ht formation from H^ and H 2 , and rates for H^ interaction with CO and 
N 2 , are taken from lPagani et al.l (|2009l) . Furthermore, all reactions involving proton transfer from 



We do not consider the neutral-neutral isotope exchange reactions proposed bv lRodgers & Charnley (2008b) 
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Hg to neutral molecules have been assigned branching ratios according to lOkal (120041) and no spin 
conversion is assumed to take place in charge exchange reactions involving . Apart from these 
special cases, reactions are assumed to produce ortho and para versions of H 2 , , and with 
equal probability. 



2.2. Other reactions 



Since nitrogen fractionation is initiated by He + reactions (see Fig. Q]), the He + chemistry 
dictates the time-scale of fractionation. In particular, the rate for the reaction 



He- 



ll, 



H" 



H + He 



(4) 



-15 cm 3 s - 



has previously been assigned widely different value s at 10 K, ranging fro m 1.1x10" 
dRodgers & Chamlevll2008bl to 2.5 x 10~ 13 cm 3 s" 1 dwalmslev et al]l2004h . We have adopted the 
measured low temperature rate Ag]=3.0x 10~ 14 cm 3 s -1 (|Schauer et al.ll 1989b . 

In order to qualitatively compare the temporal patterns in 15 N fractionation to the expected 
levels of deuteration during core evolution, we also include the chemistry for ortho and para H 2 D + 
in the model since o— H 2 can overcome the barrier for the reverse of reaction (1), i.e. 



HoD" 



o-Hc 



H^ 



HD 



(5) 



The most im portant reactions for sp in-state dependent o/p-H 2 D + formation and destruction are 
adapted from IWalmsley et al.l (|2004b . and we assume a deuterium abundance corresponding to 
HD= 1.5 x 10~ 5 relative to elemental hydrogen. 



3. Results 



Fig. |2] shows the evolution of the major chemical species. The ammonia abundance builds up 
to about 10~ 7 with respect to H 2 around a core a ge of 10 5 years (left panel), which is comparable 
to what is typically observed (IHotzel et al.ll2004h. Similar ly, the peak HCN and HNC agree well 
with values observed in dense cores (|Padovani et al .11201 1|) . At this time the CO is depleted by a 
factor of > 10 3 . A steady-st ate H 2 OPR of 2.3 x 10 ~ 5 , consistent with attempts to quantify the OPR 
in the cold dark cloud B68 (|Troscompt et al.ll2009b . is reached after about 10 6 years. 



The right panel of Fig. [2] shows the evolution of th e molecular 14 N/ 15 N r atios in the current 
model, compared to the corresponding results from the iRodgers & Charnlevl (|2008ab model. As 
expected (see e.g. Fig.Q]), the inclusion of a spin-state dependence does not significantly affect the 



-6- 



evolution of fractionation in the nitriles: the 15 N enhancement starts to build up from the initial 
pool of atomic 15 N early on, and continues to increase until gas-phase abundances have dropped to 
below 1(T 13 . The 14 N/ 15 N ratios in gas phase HCN and HNC lie in the range « 50 - 260. When 
all nitriles are frozen out, the 15 N enhancement in the bulk of nitrile ices (CN, HCN, HNC, H 3 CN) 
is ~5. 

The inclusion of spin-state dependence does however have a dramatic effect on the 15 N frac- 
tionation in ammonia. The relative delay in the timi ng of the late enhancement peak in the current 
model, as compared to lRodgers & Charnleyl(|2008ah . stems from the higher rate (~ 3 x ) of reaction 
© (see Sect. 12.21) . However, the main difference is that the NH 3 enrichment is not monotonic in 
time. The sharp decrease starting around 2x 10 5 years (time B) (i.e. an increase in 14 NH 3 / 15 NH 3 ), 
is caused by the inefficiency of reaction © as the o-H 2 abundance drops (see Fig. [2]). 

To understand what causes the fractionation patterns, we compare abundances to isotopologue 
ratios at crucial times in the evolution, marked by vertical lines A, B, and C in Fig. [2] Initially, 
the limiting reaction for fractionation is the production of 15 N + , and the 14 N/ 15 N ratio in ammonia 
is simply double that of 14 N 2 / 14 N 15 N. As gas-phase CO gets depleted after ~10 4 years (A), more 
He + becomes available to produce 15 N + , which starts to enrich the nitrogen hydrides beyond that 
of N 2 . However, after ~10 5 years, the H 2 OPR drops below 2 x 10~ 7 (B), reaction © becomes 
less efficient, and more 15 N + is instead circulated back into molecular nitrogen by 



15 M + 



N + + 14 N 2 ^ 14 N H 



+ 14 N 15 N 



(6) 



(|Terzieva & HerbstH2000T ). lowering the 15 N fraction in ammonia to less than half the elemental 
fraction. The reason for ammonia to subsequently become 15 N enriched after 10 6 years (C) is 
that the molecular nitrogen fraction decrease substantially. This both reduces the fraction of 15 N + 
converted back into molecular form (Fig.[TJ), and decreases the amount of 14 N + available to produce 
14 NH 3 . 

Fractionation changes in the bulk of amine ices follows the pattern of gas-phase fractionation, 
but with a time lag and not as pronounced extreme values. While the bulk of amine ices on grains 
has 15 N fractions within 30% of the elemental frac tion, at any given time the f ractionation in the 
outer monolayer may show much larger variations (|Rodgers & Charnlevll2008ah . 



In Fig. [3] the calculated deuterium fractionation in H 3 is plotted together with the molecular 
15 N fractionation. The D/H ratios in other gas phase molecules, in atoms, and therefore also in the 
accreted molecular ices, will closely follow that of H^ . The increase in deuteration follows very 
closely the drop in o-H 2 abundance, and after 2x 10 5 years (time B), ratios of D/H>1% are pre- 
dicted. These trends in both deuteration and OPR evolution agree well with previous investigations 
dPagani et alJboill : Slower et albood) . 



Fig. [3] shows that gas phase chemistry can indeed reproduce the requisite range of D and 15 N 
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enrichments in nitrile and amine functional groups. Three fractionation patterns develop as the 
core evolves: 



(I) For times less than about ~ 2 x 10 5 years, gaseous nitriles can develop large enrichments in 
15 N (HC 14 N/HC 15 N rj 5 - 260) but will only be modestly enriched in D, as observed (e.g. 
DNC/HNC « 0.01 - 0.05: lHirota et aU200lh : this will be reflected in the accreted ices. Ammonia 
will be similarly enriched in D but less enriched in 15 N than the nitriles. At this stage the nitrile 
gas-grain fractionation is complete. 

(II) Between ~ 2 x 10 5 — 10 6 years, the remaining ammonia can become highly fractionated in D 
but is either not enriched, or even depleted, in 15 N. 

(III) Beyond 10 6 years, ammonia becomes very highly enriched in both D and 15 N. 



Thus, the model naturally produces a range of D- 15 N fractionation in nitrile and amine func- 
tional groups that is qualitatively consistent with that found in primitive organic matter. In general, 
the highest 15 N enrichments and lowest D enrichments should be carried by nitriles. Amines 
could account for the fractionation of material that is highly enriched in D only, as well as the 
material that is simultaneously highly enriched in both D and 15 N. This functional-group de- 
pendency has already been noted in analyses of meteoritical samples. However, an important 
caveat is that conversion of nitriles to amines by H and D additio ns on grain surfaces (e.g. hydro- 
genation of HCN to CH 2 NH and CH,NHo. Eheule et alJlioill) has to be negligible. The 30% 
enrichment of the bulk ammonia ice cannot account for all the enhanced ammonia fractiona- 



tion measured in carbonaceous chondrites ( NH 3 / NH 3 ~ 185 — 256: iPizzarello & Williams 



2012[) . However, we predict 14 NH 3 / 15 NH 3 « 100 as the highest gas phase enrichment and 
this will lead to higher 14 N H 3 / 15 NH 3 ratios in the uppermost monolayers of the stratified ices 
kodgers & Charnlevll2~0~08al) . 



Interstellar gas phase 14 N / 15 N ratios in nitriles have been measured in sever al Galactic sources 
select ed for being unlikely to have experienced significant isotopic fractionation (jAdande & Ziurys 
20121) . For cold dark cores, fewer measurements exist; our model predicts very low 14 N / 15 N ratios 
for HC N and HNC. Fo r ammonia, the observed 14 NH 3 / 15 NH 3 ratios in dense protostellar cores 
340. Lis et al.ll2010f) are consistent with the modest enri chments which ca n occur within 10 5 
years. The measurement of 14 NH 3 / 15 NH 3 > 700 in L1544 (|Gerin et al.ll2009|) is consistent with 
the prediction that 15 N nuclei can be underabundant in ammonia. It also suggests a chemical age 
of ~ 3 x 10 5 years for the fractionated ammonia gas (Fig. [3 ), and that it should exist in a region 
of the LI 544 core spatially distinct from HCN and HNC (cf. IPadovani et alj|201 lb . 
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4. Conclusions 

By considering the ortho-para dependence in ion-molecule reactions involving H 2 , ammonia 
and related molecules are expected to exhibit a wide range of fractionation in both 15 N and D, 
consistent with those found in primitive Solar System organic matter. These include very large 
enrichments in both D and 15 N, as found in previous models, but also low 15 N enrichment or even 
depletion, coupled with modest to large D/H values. Nitrile functional groups on the other hand, 
while being the most likely source of 15 N hotspots, are not predicted to correlate with extreme 
D enrichment. These calculations support the premise that the 15 N isotopic anomalies found in 
meteorites were set in functional groups formed in a cold molecular cloud. 

The predicted 14 N/ 15 N and D/H ratios can account for the limited measurements existing 
for interstellar cloud cores. To further test and constrain these theoretical models, surveys for 
key molecules and measurements of relevant 14 N/ 15 N ratios in molecular clouds are required. 
A more extensive analysis of how correlations between deuterium and 15 N enrichment in specific 
molecules are affected by spin-state-dependent chemistry is beyond the scope of this Letter and will 
be considered elsewhere (Wirstrom et al. 2012, in preparation). Predictions from these models can 
be tested by mea surements of various multiply-substituted isotopologues (e.g. DC 15 N, 15 NH 2 D, 



Gerin et al.ll2009l) . which will be facilitated with the Atacama Large Millimeter Array (ALMA). 



This work was supported by NASA's Origins of Solar Systems Program and the Goddard 
Center for Astrobiology. 
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Fig. 1. — Chemical network showing the main reactions responsible for 5 N enhancement in ni- 
triles and ammonia. 
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Fig. 2. — Left panel: Time evolution of the nitrogen chemistry in dense cores, compared to CO and 
o-H 2 . Crucial times for the 15 N fractionation in nitrogen hydrides are marked by vertical dotted 
lines, labelled A, B, and C. Right panel: Evolution of the 14 N/ 15 N ratio in gas-phase molecules 
and bulk of ammonia and nitrile ices. Black curves r epresent the current model, w hile red curves 
show the corresponding isotopic ratios in the model of lRodgers & Charnleyl (|2008af) . The assumed 
elemental value of 440 is marked by a dotted horizontal line. Values for gas-phase nitriles are 
omitted after ~2xl0 5 yrs when abundances for both isotopologues are essentially zero due to 
freeze-out. Dashed curves represent ice-phase species. 
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Fig. 3. — Nitrogen fractionation in ammonia and HCN, both gas-phase (solid) and bulk ice 
(dashed) in blue curves according to left-hand axis, compared to deuterium fractionation in 
H 2 D + /Hg (green curve and scale on right hand axis). 15 N enhancements are given relative to the 
elemental 15 N/ 14 N ratio of 1/440, while D enhancements are relative to elemental D/H=l .5 x 10~ 5 . 



